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We introduce and numerically solve a mathematical model of the experimentally established 
mechanisms responsible for the symmetry breaking transition observed in the chiral crystallization 
experiments reported by I. Weissbuch, L. Addadi, L. Leiserowitz and M. Lahav, J. Am. Chem. 
Soc. 110 (1988), 561-567. The mathematical model is based on five basic processes: (1) The 
formation of achiral glycine clusters in solution, (2) The nucleation of oriented glycine crystals at 
the air/water interface in the presence of hydrophobic amino acids, (3) A kinetic orienting effect 
which inhibits crystal growth, (4) The enantioselective occlusion of the amino acids from solution, 
and (5) The growth of oriented host glycine crystals at the interface. We translate these processes 
into differential rate equations. We first study the model with the orienting process (2) without (3) 
and then combine both allowing us to make detailed comparisons of both orienting effects which 
actually act in unison in the experiment. Numerical results indicate that the model can yield a 
high percentage orientation of the mixed crystals at the interface and the consequent resolution 
of the initially racemic mixture of amino acids in solution. The model thus leads to separation 
of enantiomeric territories, the generation and amplification of optical activity by enantioselective 
occlusion of chiral additives through chiral surfaces of glycine crystals. 



I. INTRODUCTION 

Theoretical proposals for prebiotic chemistry suggest that homochirality emerged in nature in abiotic times via 
deterministic or chance mechanisms [l| . The abiotic scenario for the emergence of single homochirality in the biological 
world implies that single asymmetry could have emerged provided a small fluctuation from the racemic state can be 
amplified to a state useful for biotic evolution. For this reason, experiments that can demonstrate the feasibility of 
stochastic mirror symmetry breaking involving the self assembly of molecular clusters, and in possible conjunction with 
interface effects, are particularly important. This is because, once generated by chance or an initial chiral fluctuation, 
the chirality can then be preserved and transmitted to the rest of the system provided that the symmetry breaking 
step is coupled to a sequential step of efficient amplification via self-replication reactions. Some relevant features 
common to such systems are that they take into account the small fluctuations about the racemic state and that they 
display nonlinear kinetic effects. Stochastic scenarios are theoretically well understood^, 0] and are experimentally 
feasible in the laboratory |4[ . 

Among the experiments dedicated to exploring chance mechanisms in chirality, a particularly noteworthy and 
important result with a marked relevance for prebiotic chemistry stands out. The experiment we are interested in 
modeling here was reported some years ago by the Rehovot group which dealt with an autocatalytic process for the 
resolution of racemic a-amino acids within crystals of glycine grown at the air/solution interface They applied 

cooperative crystallization processes for the spontaneous separation of racemic mixtures of a-amino acids rich with 
glycine into optically pure enantiomers. Their experimental model involves slow evaporation of aqueous solutions of the 
ccntrosymmetric form of glycine containing racemic mixtures of a-amino acids. Due to the unique crystal structure of 
glycine, all chiral D-amino acids except for proline are occluded within the crystal through the (010) face, whereas the 
L-amino acids are occluded through the (010) face. Glycine crystals float in solution, so that only one face is available 
for growth. Thus when glycine crystals are grown at the air/water interface in the presence of DL-amino acids, only 
one of its enantiotopic faces, say (010), is exposed to the solution and so picks up only the D-amino acid together 
with glycine. By symmetry, crystals exposing their (010) face towards solution occlude only the L-enantiomers. Now, 
if by chance a single or small number of oriented glycine crystals grow initially at the interface, the bulk solution will 
be enhanced with the amino acid of one handedness. The preservation and transmission of the chirality generated by 
chance of the original Adam 23\ crystal means that new crystals grown at later stages at this interface must adopt 
the same orientation. There are two proven ways this is achieved^, @: by means of an (i) hydrophobic effect and (ii) 
by a kinetic inhibition effect. Regarding the first effect, if the solution contains hydrophobic amino acids, these tend 
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to accumulate at the interface forming two-dimensional domains acting as templates for the oriented crystallization 
of the glycine crystals. Thus, the L-amino acids induce crystallization of floating glycine crystals exposing their (010) 
face towards solution and these occlude only the D-amino acids. This asymmetric induction has been established 
experimentally^, 0|- As for the second effect, this comprises an enantioselective inhibition of the glycine nuclei by 
the amino acids present in solution (these can be both partially dissolved hydrophobic as well as hydrophilic amino 
acids). This independent effect was proven experimentally by achieving complete orientation of the floating glycine 
nuclei when grown in the presence of DL-leucine and hydrophilic L amino acids. The presence of the DL-leucine is to 
ensure nucleation of floating glycine crystals exposing either enantiotopic face towards solution in a 1 : 1 ratio whereas 
increasing the concentration of the hydrophilic amino acid additives inhibited the glycine nuclei exposing their (010) 
face and so prevented their further growth. These hydrophilic amino acids inhibit the crystal nuclei from growing 
and developing into macroscopic crystals Q. Both the hydrophobic and kinetic effects act cooperatively in the same 
direction in that they both contribute to the territorial segregation of the enantiomers. 

The overall experimental process can be summarized by the following steps (Figure [!} [6] , 4|:(i) mirror symmetry 
breaking via the enantioselective occlusion of one of the enantiomers of the racemic a-amino acids within crystals 
of glycine grown at the air/solution interface, (ii) the oriented crystal formed at the interface operates as a seed for 
further occlusion of amino acids of the same handedness, (iii) the amplification experiments comprises self- aggregation 
of hydrophobic a-amino acids into chiral clusters that operate as templates for an oriented crystallization of fresh 
crystals of glycine, (i.e, the hydrophobic effect), and (iv) enantioselective inhibition of embryonic nuclei of glycine 
generated at the air/solution interface by the water soluble a-amino acids formed in excess during the process, (i.e., 
the kinetic effect). The fundamental relevance of this experiment is that it provides a simple chemical model for the 
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FIG. 1: The overall process of mirror symmetry breaking and amplification. Scheme adapted from references [J, [6j]. 

generation and amplification of optically active amino acids in prebiotic conditions. They used a centrosymmetric 
crystal structure of glycine, which is the only achiral amino acid and is the major component found in modern prebiotic 
synthesis experiments The possibility to maintain and propagate the same chirality to the entire system can lead 
in principle to a complete separation of enantiomeric territories and subsequent optical activity. The enantiomeric 
resolution and symmetry breaking is achieved without the need to input mechanical energy: in particular there is 
no need for stirring [10j nor grinding [111 ], and it involves only amino acids which have an immediate and obvious 
relevance for prebiotic chemistry. That the experiment dispenses with the need for mechanical energy is a most 
remarkable feature. Some provisional explanations of how this is achieved are offered in the final conclusions. 

In spite of its importance, we are unaware of any prior attempts to offer a mathematical model for describing the 
reported phenomenon. The purpose of this paper is to provide a simple and minimal mathematical model of the key 
processes to further elucidate the mechanisms responsible for the observed symmetry breaking and to gain further 
understanding of the experiment itself. Our aim is to describe the resolution process by means of a kinetic mean field 
model in which we include only the minimal essentials responsible for the crystal growth and the symmetry breaking 
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phenomenon. The mathematical model is based on five basic processes that parallel closely those that have been 
established experimentally: (1) The formation of achiral glycine clusters in solution, (2) The nucleation of oriented 
glycine crystals at the air/water interface in the presence of hydrophobic amino acids, (3) A kinetic orienting effect 
which inhibits crystal growth from embryonic crystals, (4) the enantioselective occlusion of the amino acids from 
solution, and (5) the growth of oriented host glycine crystals at the interface via glycine up-take. We first introduce 
the model with only the hydrophobic orienting process (2) as an approximation to the complete model. This is a 
justifiable first approximation to consider given that the experimental results prove that this effect plays a dominant 
role in the orientation of the glycine crystals |6J. The major simplifications are introduced at this stage and we use 
this model to prove that the symmetric state is unstable. We then include both effects (2) and (3) in a complete final 
version allowing us to make detailed comparisons of both effects which actually act in unison in the experiment. We 
present numerical studies of the effects of varying amounts of both the hydrophobic and hydrophilic additives on the 
orientation of the floating glycine crystals, as was originally considered in the experiments. The results indicate that 
our simple model can yield a high percentage orientation of the glycine-host plus guest-enantiomer mixed crystals 
at the interface with the consequent resolution of the initially racemic mixture of amino acids in solution. The 
mathematical model thus leads to separation of enantiomeric territories, the generation and amplification of optical 
activity by enantioselective occlusion of chiral additives through chiral surfaces of glycine crystals. 



A full experimental account is given in the two papers [a, |6| and is reviewed[4j, [l2j together the references therein; 
and these serve as basic inspiration and background guide for elaborating the model we present in this paper. 

Following the experimental conditions, we consider a supersaturated solution of glycine and racemic a amino acids. 
The glycine monomers and clusters in solution will be denoted by A\ and A r , respectively, the cluster being made 
up from r-monomer units. It is reasonable to assume that a typical glycine cluster size should exist < r >= M for 
M monomeric units. This is the average size of the primary nucleation or seed crystal which can then precipitate at 
the air/water interface where it will, in the presence of hydrophobic D and L-amino acids, assume one of two unique 
orientations jgj. Let Xm and Ym denote the two possible orientations of the glycine crystal floating at the interface: 
X stands for the (010) enantiotopic face exposed to solution and Y for the other (010) enantiotopic face exposed to 
solution. The D and L denote the two amino acid enantiomers. For a strictly racemic solution of the amino acids, 
(which is of course impossible to realize experimentally 1 1 3j ) . we would expect an equal X : Y = 1 : 1 proportion 
of the two glycine crystal orientations at the interface. However, if for example there is an imbalance in the initial 
concentrations, say L > D, then more of the F-oriented crystals will be precipitated. Likewise for the X-oriented 
crystals if instead L < D. 

Next, we consider the enantioselective occlusion of the amino acid into the host glycine crystals: recall only the 
(010) face of the glycine crystal can occlude the D amino acid whereas only the (010) face can occlude the L amino 
acid; see Figure [T] It will prove instructive to first write down a kind of "microscopic" model which contemplates 
all kinds of crystal aggregations growing from both glycine take-up and the amino acid occlusion and with cluster 
size dependent rates and then simplify by taking the reaction rates independent of cluster size, and assuming only 
cluster-monomer aggregation. We will not consider fragmentation in any case: we assume irreversible steps from the 
outset. Since the hydrophobic effect dominates over the kinetic one@, we will consider this one first. 



We now introduce the explicit processes to be included in the mathematical model. We transcribe the scheme in 
Figure [1] into reaction steps. Simplifications will follow later on in order to obtain what we will consider as a minimal 
model leading to mirror symmetry breaking. 

The formation of glycine aggregates/clusters in solution. The achiral molecules of glycine in solution combine 
pairwise to yield achiral glycine clusters in solution, where r, s denote number of monomers in the cluster and S r , s is 
the cluster size dependent rate constant: 



II. REACTION MODEL 



A. 



Preliminaries 



B. Model based on the hydrophobic crystal orientation effect 




(1) 



We assume this step is irreversible, there is no fragmentation of the glycine aggregate. Below, we simplify this so that 
only dimers are formed (r = 1, s = 
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Nucleation of oriented glycine crystals at the air /water interface in presence of hydrophobic LD -amino acids. The 
hydrophobic amino acid at the interface acts as a nucleator or seed (see also Scheme 5 in [J]): 

A r + D 1 {DiY r }. (2) 

Here and henceforth, the brackets denote the solid phase. This process occurs at the rate fi r . The two enantiotopic 
crystal faces exposed to the solution are denoted by X r and Y r , respectively. Realistically, a single leucine or valine 
molecule cannot serve as a nucleus for the glycine crystallization, instead, several dozens arranged in self-assembly 
are usually required. However, for the sake of simplicity (and limitations on computational time), we assume that one 
hydrophobic amino acid monomer is sufficient to trigger the nucleation. Otherwise, we would have to additio nally 
model the formation of the chiral template clusters formed by several hydrophobic amino acids at the interface 0, Il5j . 
This template, whatever its actual size, is incorporated into the growing glycine crystal, but in the face exposed to 
air. We note that HPLC analysis of the crystals indicates that only minute amounts of the orienting hydrophobic 
aminos acids are occluded through the upward pointing faceQ, so our single monomer "template" assumption is not 
unreasonable and results in a welcome mathematical simplification. Varying the size of the template can only have 
quantitative but not qualitative effects in so far as mirror symmetry breaking is concerned. Below, we will assume a 
transition regime, that is, only dimcr glycine clusters r = 2 get nucleated in this way to form oriented host seeds. 

Enantioselective occlusion of the amino acid monomers from solution. Incorporation of the guest molecules (amino 
acids) into mixed host crystal leading to formation of homochiral mixed crystal. As dictated by the actual experimental 
results Q, we assume that the X face exposed to solution occludes only the D-amino acids, whereas the Y face occludes 
only the L-amino acids: 

{L x X r D n } + D x ^ {L x X r D n+1 }, 

{D{Y r L n } + L\ H {D x Y r L n+x }. (3) 

The corresponding rate £ rj „ could depend on both the size r of the glycine host as well as on the number n of 
previously occluded enantiomers. We assume the oriented host glycine crystal occludes one amino-acid monomer 
at a time. Crystallographic models suggest that single molecules, not clusters, are incorporated one at a time into 
the growing crystal. The notation in {...} is such that reading from left to right: the orienting hydrophobic amino 
acid attached to the face exposed to air, the enantiotopic crystal face exposed to solution (composed of r-glycine 
monomers) and the number n of occluded amino acids from solution. There is no structural or sequence information 
implied. Note this step implies that the percent ee of occluded enantiomers per host crystal will be 100%, since each 
enantiotopic crystal face occludes enantioselectively. This is fully justified by the actual experimental results, see the 
fourth column of Table I in reference @. 

Growth of the oriented host glycine crystal. Growth of the enantiotopic face exposed to solution by take-up of the 
achiral glycine monomers from the solution: 

{LxXrDn + A X ^ {L X X r+X D n }, 

{D x Y r L n } + Ax ^ {D x Y r+x L n }. (4) 

The newly acquired monomers adopt the same orientation as their host-face. Assumption: take-up involves one 
glycine monomer at a time. Here the rate a r ^ n could in principle depend on r and n. 

The above steps give a specific articulation of the Scheme represented in Figure [TJ The sequence of pictures in 
Figure [2] illustrates how the the hydrophobic orienting effect plus the enantioselective occlusion work in tandem to 
yield resolution and optical activity. For illustrative purposes only, we consider a model in which one hydrophobic 
amino acid is sufficient to orient crystals at the air/water interface (see the above remarks) and that the crystals can 
occlude up to two amino acids from the solution. The events are ordered from left to right as indicated by the arrows. 
Start from a closed system containing a pool of amino acid enantiomers in an initial ratio of L : D = 10 : 9 plus an 
initial concentration of glycine aggregates, Gly. By chance, a glycine aggregate nucleates with its (010) exposed to 
the solution, this step uses up one hydrophobic L from the solution so now we have L : D = 9 : 9. The crystal starts 
occluding the D monomers from a solution and it therefore becomes enriched in the L enantiomer: L : D = 9 : 7. It is 
more probable that the next crystal to be nucleated at the interface will have the same orientation as the Adam crystal. 
Another L is needed to orient the crystal and so two more D monomers get occluded. The solution is further enriched 
in the L enantiomer. In the course of time, this yields a cascade mechanism and results in a final configuration with 
a fully oriented glycine crust with a resolved amino acids for each individual crystal and an optically active solution 
rich in the L-enantiomer. The finite surface area gets covered by nucleated glycine crystals, and the reactions at the 
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interface cease when either the maximum number of -D's per crystal are occluded or when they are depicted from the 
solution. 

Note in passing: there no explicit spatial or coordinate dependence in this model. We do not explicitly distinguish 
the interface surface from the bulk solution. We can think of this as a homogeneous two-dimensional model, namely, 
the two-dimensional air/water interface, with the role of the three-dimensional bulk solution (below the interface) as 
providing the interface with the new glycine and amino acid molecular building blocks. Also, there is no fragmentation; 
once formed the crystals do not break up into smaller pieces, the guest molecules stay with their hosts. The real 
experimental situation is three-dimensional, being composed of an interface plus bulk solution. It is clear that the 
molecules in solution must diffuse in order to reach the interface. Even the crust formed at the interface is not 
homogeneous, but is composed of the fusion of many plate- like crystals @. We ignore size distribution of the crystals 
grown at the interface. But for the purposes of obtaining a simple model that leads to experimentally established 
resolution, these details will not matter. When we treat the complete truncated model below, we will however 
implement a way to effectively account for the finite surface area available for the oriented glycine crystals. 
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FIG. 2: Simplified schematic representation of our model for the generation and amplification of optical activity by enantiose- 
lective occlusion of amino acids by floating glycine crystals oriented at the air/water interface by hydrophobic additives. The 
system starts from an initial enantiomeric excess L > D in solution. The horizontal arrows indicate the temporal sequence of 
events. 



III. MATHEMATICAL MODEL 



We transcribe the above reaction processes into the corresponding differential rate equations. Introduce concentra- 
tion variables a r (t) = [A r ], L\(t) — Di(t) = [Di] and for the two host-plus-guest crystal orientations exposed 
towards the solution f r ,n{t) = {LiX r D n } and f r , n (t) = {DiY r L n }, respectively. We consider separately the kinetic 
equations for f r fi, /2,n, and then / r>n and similarly for the other enantiotopic face orientation. That is, we can treat 
the growth of pure host crystal, growth of guest on minimal host substrate and then the general mixed growth of both 
host and guest. Remember f r ^ — f host, guest, so there must be a minimum host size r to accommodate the occluded 
guests. 



6 



Applying law of mass action we obtain the following kinetic equations: 

oo 

-fiia^Lt + Di) - 28na\ - ^ 5i.ka\a k 



dai 
~~dt 



fe=2 



X! X! Olr,n{fr,n + fr,n), (5) 



r n=0 

r-1 



da r 1 ^-^ 

fe=l 



d£>i 

dfrfi 

dt 

dfrfi 

dt 
dt 

dfr, n 

dt 



-^25 r%k a r a k , (r > 2) (6) 
fe=i 

Mr a r ~ Li £,r,nfr,n, (?) 

r r n— 

-A. ^ /X r a r - -Dl ^ ^2 £r,nfr,n, (8) 
r r n— 

+ai(a r ._i,o/r-i,o - a r;0 /r,o), (9) 

H r D\a r — irflfrfiLi 

+a\(a r -\fif r -ifl - a r ,o/ r ,o), (10) 

Dl(,^r,n—lfr,n—l ^r,nfr,n) 

+ a 1 (a: r -l,n,/r-l,n - a r , n /r,n), (« > 1) (H) 

^1 (£r,n— 1 fr,n— 1 ^r,nfr,n) 
+ai(a r _i, n / r _i /r,n), (« > 1). 

(12) 



For a closed system, we have the constant density constraint: 

oo oo oo 

Li+L>i+5Z ra '' + Z!IZ( r + n + 1 )(^« + ^«) =consi ' ( 13 ) 

r r n 

The main problem with such a model is the large number of parameters that have been introduced: namely 8 r ^ ni fi ri o^ r n 
and £ r , n . It is thus clear we need to make several simplifications already at this stage. 

A. Simplifications of the model 

We assume that only the coalescence of two glycine monomers A\ to form an A2 dimer needs to be retained, and 
only A2 dimer clusters are used to form new oriented host crystals when they encounter a hydrophobic amino acid 
monomer. We also assume that all rates are independent of cluster size, so that 

#1,1 = 5, 5 r ,k — 0, {otherwise) (14) 

[12 = (J., [ir = 0,(r = lkr > 3) (15) 

a r , n = a,(r > 2), ai ; „ = (16) 

£r,n = £,(r>2), £l,n = 0. (17) 



Then the above kinetic equations take the following form: 

-25a\ - am £(/r,n + fr,n), (18) 
r=2 n=0 

5a\- /^(Li+Di), (19) 

-nLia 2 - £L-l y~] fr tn , (20) 

r=2 n=0 

-/iDia 2 -^Di25Z/ r ,„, (21) 

r=2 n=0 

HL 1 a 2 -(£D 1 +aa 1 )f 2 fl, (22) 

Ml>iO2-(Cii+aoi)/2,0, (23) 

-£Di/r,o + cmi(/r-i,o - /r,o), (r > 3) (24) 

-£&i/r,o + aai(/ r _i, - /r,o), (r > 3) (25) 

^l(/2,n-l - /2,n) - aai/ 2 ,„, (« > 1) (26) 

£M/2,n-i ~ /2,n) - aai/ 2 , n , (n > 1) (27) 

£Dl(/r,n-l - /r,n) + ™l(/r-l,n - /r,n)j 

(r>3,n>l) (28) 

£-ta(/i-,ra-l — /r,n) + £*Ol(/r-l,n — /r,n)j 

(r>3,ra>l). (29) 



da! 

~dt~ 

r=2 n=0 

d02 _ x „ 2 

dt 

dDi 
dt 

df 2,0 

dt 

df 2,0 

dt 

dfrfl 

dt 

dfrfl 

dt 

dj%n 

dt 
dt 

df r,n 

dt 

dfr.n 

dt 



As an important check on internal consistency, we verify explicitly that mass conservation holds. To this end, we 
consider the following change of variables 

oo oo oo oo 

P = frfi, M 2 = M = J2J2 (30) 

r — 3 n—1 r — 3 n—1 

oo oo oo oo 

P = frfi, M 2 = Y. hn, M = Yl E fr,n- (31) 



r— 3 n—1 r— 3 n— 1 
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Then the above infinite collection of equations can be written in terms a mathematically closed set of just twelve 
equations: 

da i 



dt 



-25a\ 

-crai (/ 2 ,o + M 2 + P + M + f 2fi + P + M + M 2 ) , 





(32) 


5a\ - ^ a 2{L\ + Di), 


(33) 


-pL x a 2 - £Li (/ 2)0 + P + M + M 2 ) , 


(34) 


-liD x a 2 - iD x (/ 2 , + M 2 +P + M) , 


(35) 


pLia 2 - (££>i + aai)f 2 . , 


(36) 


pD x a 2 - (£Li + aai)/ 2 ,o, 


(37) 


-£DxP + aai/2,0, 


(38) 


-^Z/i-P + aai/2,0, 


(39) 




(40) 


££1/2,0 - aaiM 2 , 


(41) 


^DiP + aaiM 2 , 


(42) 


^iiP + aaiM 2 . 


(43) 



da 2 

~dT 

dDi 

<#2,0 

dt 
df 2 ,o 
dt 
dP 

~dt_ 
dP 

~dt 

dM 2 

dt 
dM 2 

dt 
dM 

dt 
dM 
~dt 

We need to have kinetic equations for the corresponding densities. The following definitions will accomplish this: 



p P = X> + l)/r,o, (44) 

r=3 
00 

PM2 = 5Z(n + 3)/ 2 , n , (45) 

71=1 
OO OO 

pM = ^^(r + n+l)/ r>n . (46) 

r=3ra=l 

Then the kinetic equations for the densities follow: 

PP = -^D 1 p P + aa 1 (4/ 2 ,o + P), (47) 
PM2 = 4££>i/ 2 , + ££>iM 2 - aaip M 2, (48) 

Pm = ££>i (-P + Pp + A'/) + aoi (p M2 + M 2 + M ) . (49) 

There are of course analogous equations for the Z 2 or chiral partners p~p,pM2 and p~M- To express the mass conservation 
in terms of these densities, define 

P = 3/2,o + PM2 + Pp + Pm, (50) 
P = 3/2,0 + PM2 + Pp + Pm, (51) 

then the following constraint holds 

ai + L\ + D\ + 2a 2 + p + p = const. (52) 
This proves that total mass of the reactants is conserved in our scheme. 
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B. Linear Stability 

It is useful at this stage to check if the model is actually capable of mirror symmetry breaking. Consider the 
following variables where H denotes pure host crystals, and M the mixed host plus guest crystals: 

H = f 2 ,o + P, (53) 

H = hfi + P, (54) 

M = M 2 + M, (55) 

M = M 2 + M. (56) 

Then we can further reduce the system from twelve to only eight differential equations: 

da i 

~~dt 

da 2 

~~d~t 
&L X 

~dT 
dDi 

~~dT 
dH_ 

~~dt 
dH 

~~dt 

dM 

dt_ 
dM 
~dT 

We can define three types of chiral polarizations ry, 9, <fr. Investigate mirror symmetry breaking by the change of 
variables: 

X = Li+Dx,w = H + H,z = M + M, (65) 

Li-Di a H — H M-M 
V = ,8= A= ■ (66) 

X w z 



-28a\ -aax(H + H+ M + M) , 


(57) 


Sal ~ ^ a 2{L\ + Di), 


(58) 


-hL^-zlJh + m), 


(59) 




(60) 


fiLia 2 - £,D 1 H, 


(61) 


liD\a 2 — ^L\H, 


(62) 




(63) 




(64) 



The transformed equations Eqs. ([57H6"4"]) read as follows: 

-28a\ — aa x (w + z), (67) 

Sal ~ ^ a 2X, (68) 

'^a 2 x - |x(w(l - rj6) + z(l - rj<f>)j , (69) 



da\ 
~~dt 

dt 

dt 

§ = f(l-7? 2 )H + ^), (70) 

— = /ia 2 x —(l-r)6), (71) 

d0 ^a 2 x ( „, £,X n a2 , /vo x 
dt w 2 

dz _ £xw 

dt ~ 2 



(1 - rfi), (73) 



"4 = ^(0-v)-^d-v0). (74) 
dt 2z zz 
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Consider the stability of the symmetric solution r\ = 9 = (f> — 0. Linearize the system of the three corresponding rate 
equations and determine stability of the solution from the Jacobian matrix: 



Expansion along the minors yields the determinant det which is also given by the product of the three eigenvalues Ai : 



This result immediately tells us that the symmetric or racemic state is unstable. To confirm this, first consider the 
case where all three eigenvalues A^ are real. Then the only way for their product to be positive is if the individual 
eigenvalues have the algebraic signs (+,+,+), or (+,—,—) (the order is not important): the point is, either one 
eigenvalue or else all three must be positive, and thus in both cases the symmetric state is unstable. So chiral 
symmetry will be broken in this infinite but closed model. Note that the remaining sign possibilities (—,—,—) and 
(+, +, — ) (again, the order is not important) would yield instead a negative determinant det < 0, and could imply 
stable or unstable, respectively. This latter situation is inconclusive without further analysis. If on the other hand 
det has complex eigenvalues, these always occur in complex conjugate pairs, and without loss of generality we may 
suppose that A3 = A2*. Then det = Ai | A2 1 2 > if and only if Ai > 0, and again we conclude that the symmetric state 
is unstable. The expression for det depends on the rate of enantioselective occlusion £ and the rate for hydrophobic 
orientation of host crystals /i and suggests that the former is more important than the latter for driving the symmetry 
breaking: this is so because we can set /x = in Eq. (|76[) and det > remains positive, so the racemic state remains 
unstable. But if we set £ = 0, then det vanishes identically and the result is inconclusive at this lowest order. On 
the other hand, we note that neither the rate of glycine up-take a nor the rate of glycine cluster formation 6 appear 
explicitly in the determinant, suggesting that these processes in and of themselves are not the decisive factors for 
symmetry breaking. These expectations are in accord with the processes depicted in the scheme represented in Fig. [T] 
the crucial amplification cycle involves the orientation and occlusion steps only. The processes of glycine aggregation 
in solution (5) and growth by glycine up-take (a) are prior events that are not involved in this cycle. Thus, at this 
stage, we may be confident that our model will qualitatively capture the main symmetry breaking aspects of the 
actual experiment, and that the bare minimal ingredients required for this are crystal orientation and enantioselective 
occlusion. 

The above model allows unlimited growth of the mixed host plus guest crystals. A severe truncation yielding a 
minimal model consists of an r = 2 glycine host with at most one occluded amino acid guest n — 1 per host. The 
determinant in this case is given by Eg. (1751) after deleting the t^^w/di term. The linear stability analysis carried out 
for this truncated model leads to the same conclusions as above. We thus find that the symmetric state is unstable 
in two opposite limiting cases of the underlying model: for glycine host dimers and one occluded guest as well as for 
arbitrarily large glycine host crystals r < 00 accommodating an unlimited number n < 00 of occluded guests. Thus 
we may expect that the symmetric state will be unstable for all intermediate host/guest truncations of the model. 
This expectation is confirmed by the numerical simulations presented below. 



The model leads to the unlimited growth of both the glycine host crystals as well as to an unlimited number of 
enantiomers that can be occluded per host crystal. The experiments by contrast yield finite size mixed glycine plus 
host crystals with a rather uniform size distribution and with each host occluding a small percentage of the available 
amino acids [HQ. We must also recognize the fundamental physical limitation imposed by the fact that the air/water 
interface is of finite area: the experiments are carried out in a bounded reaction domain which means only a finite 
air/water interface is available for the crucial orientation/amplification processes to take place. We will thus need to 
model the effect of a bounded (finite area) interface. Once the interface is covered by floating glycine crystals, no more 
hydrophobic amino acids can diffuse up to the surface, the remaining reactions from this point on can only be the 
glycine up-take from the solution and the enantioselective occlusion of the amino acids in solution. The finite size of 
the experiment implies that these latter two processes must be limited as well; we will thus consider how to truncate 
both the host growth as well as the number of occluded guests per host. The reactions will thus stop at some point: 
there must be maximum values of both R and N such that 2 < r < R and hosts < n < N(R). Each instantaneous 
size r of the host crystal will be allowed to occlude up to maximum number n(r) of guests. We can account for all 




(75) 




(76) 



i=l 



IV. COMPLETE TRUNCATED MODEL 



11 



of these limiting features by implementing certain truncations or cut-offs that we apply to the underlying model. At 
the same time, we complete the model by including the kinetic inhibition effect. 



A. Hydrophobic and kinetic effects and hydrophobic and hydrophilic additives 

The experiment demonstrates that the glycine crystals at the interface are oriented via two distinct effects: (i) 
hydrophobic and (ii) kinetic. The former is due to the induction by hydrophobic amino acids while the latter is 
achieved through the inhibition of nucleation and growth of the oriented crystals. The hydrophobic effect is confined to 
hydrophobic amino acids but the kinetic effect applies to all a-amino acid additives, both hydrophobic and hydrophilic. 
It is important to bear in mind that there is hence not a one-to-one correspondence between effects and additives. For 
this reason, in the full model below, we introduce the two kinds of additives and talk about the effect of the additives 
rather than the hydrophobic or kinetic effects per se, since the latter are not easily separable. The two effects however 
do act in the same direction. The inclusion of both types of additive leads to only a minor modification and leads to 
a complete model. We also implement a consistent truncation where R is maximum size of the glycine host crystal 
and n = n(r) is a maximum number of occluded amino acid guest monomers for any instantaneous host size r. 

For completeness, we list the full set of chemical reactions defining the final model. Formation of glycine crys- 
tals/clusters in solution. 

(77) 

Diffusion of amino acids in bulk solution to interface implying nucleation of oriented glycine crystals at the air/water 
interface in presence of hydrophobic LD-amino acids. 

A 2 + L, ^ {L^}, 

A 2 + D 1 ^ {D1Y2}. (78) 

The effective conversion rate of bulk amino acids to amino acids at the interface will depend on a critical glycine 
crystal interface concentration f c . Beyond this concentration, no more hydrophobic amino acids will diffuse up to the 
surface, and consequently no more fresh glycine crystals can nucleate at the surface. To implement this we will take 

M/c) = tt>e(/ c -/(t)), (79) 

where f(t) is be the instantaneous concentration of glycine crystals, with and without occluded amino acids, at time 
t, and f c is a critical surface concentration which is supposed to effectively mock the finite interface. Here: 

R r 

H(t) = M e(/ C - £(/r,n(*) + /r,n(*)) , (80) 
r— 3 n— 

where the unit step function is denned as 

e (*)={J;*<o° ( 8l ) 

Once the instantaneous concentration of glycine crystals reaches f c then this reaction shuts off, from which point 
on the hydrophobic orienting effect ceases to act. This is an effective way of implementing the finite interface area 
constraint without adding complicated spatial dependence to the model. 

In keeping with the above remarks and experimental facts, we must allow for the hydrophobic amino acids to inhibit 
those glycine nuclei exposing their enantiotopic faces towards the solution from further growth. We also introduce a 
second species of hydrophilic amino acids L, D, which can only participate in the kinetic effect. 

Inhibition of the fresh glycine seeds by kinetic effect due to the hydrophobic additives. 

{L 1 X 2 } + D 1 -!* {L 1 X 2 D 1 }, 

{DiY 2 } + L\ {£>iy 2 £i}. (82) 
Inhibition of the fresh glycine seeds by kinetic effect due to the hydrophilic additives: 

{L 1 X 2 } + D 1 -A {L^M, 

{D 1 Y 2 } + L 1 A {£>iF 2 ii}, (83) 
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the total inert product of this inhibited crystal seeds (via either the hydrophobic or hydrophilic additives) will be 
denoted as P(t). An excess of L > D in solution inhibits glycine nuclei exposing their (OlO) towards solution^, 
preventing their further growth, and analogously for D > L. We have opted to model this as a cross-inhibition 
reaction, in the spirit of the Frank model, except here we form a small tetramer unit with a glycine dimer trapped 
between the hydrophobic templater and the hydrophilic growth inhibitor. The templater and growth inhibitor have 
opposite handedness. So, if the glycine host crystal occludes an amino acid having the opposite chirality of the 
nucleator amino and located on the opposite enantiotopic face, then the seed glycine crystal has been inhibited and 
can not grow further. Otherwise, if the seed grows by incorporating a glycine monomer, it can escape this "bottleneck" , 
and will continue to grow via glycine up-take or by amino acid occlusion. 

The remainder of the reactions are unchanged, except that now we impose the finite size truncations (bounds on r 
and on n) remarked above. Growth of the oriented fresh glycine hosts. 

{L 1 X 2 }+A 1 {£1X3}, 

{D 1 Y 2 }+A 1 -A {D 1 Y 3 }. (84) 

Growth of the oriented host glycine crystal. 

{LiXrDn} + {L x X r+1 D n }, 
{DiY r L n } + {DiY r+1 L n }, 
(3 < r < R - 1,0 < ri < T(r)). 

(85) 

Enantioselective occlusion of the amino acid monomers from solution. 

{L x X r D n } + D X ^ {L x X r D n+1 }, 

{DiY r L n } + Li {DiY r L n+l }, 
(3 <r < R,0 <n < T(r) - 1). 

(86) 

Each host glycine crystal can grow by incorporating achiral glycine monomers from the solution, so we must impose 
a maximum number R of glycine monomers forming the host crystal. By the same token, we impose a maximum 
number of guest chiral monomers N to be occluded in the host crystal. The maximum number of guests that a crystal 
can occlude depends on the number of host monomers of the crystal, so we define T(r) = Floor[yr], where Floorer] 
rounds the elements of -fr to the nearest integer less than or equal to and < 7 < 1 is a free parameter that allows 
us to vary the percentage of occluded amino acids. 
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With these truncations and limits and for both types of additives, the final set of kinetic equations reads as follows: 

-28a\ 



da i 
~dt 



R-lT(r) 

-aOi(/ 2 ,0 + ho + S(/r,n + fr,n)), (87) 

r— 3 n— 

^ = Sa\-iMn[Lx + D x ), (88) 



eft 

eft 

~dT 
dP 
~dF 

dfr.n 

dt 



R r(r)-l 

-Miia 2 -^i(/2,o + ^ E ( 89 ) 

r — 3 n— 

fl r(r)-i 

r=3 n=0 

-PhflLu (91) 

-0h,oRi, (92) 

£(£1/2,0 + £1/2,0) + £(£1/2,0 + Li/2,0), (93) 

A(r, n)/zLia 2 + S(r, n)£-Di/ r , n _i - C(r, n)£Di/ r ,n 

+£>(r, n)aai/ r _i ; „ - £(r, n)aoi/ r ,„ - F(r, n)/3L>/r,™, 

(2<r<R,0<n<T(r)) (94) 

A(r, n)^£>ia 2 + B(r, ra)££i/r,n-i - C(r, ra)££i/r,n 

+D(r, n)aai/ r _i ; „ - £(r, n)aai/ ri „ - F(r, n)f3Lf r . n 

{2<r <R,Q<n<T{r)). (95) 

For a more streamlined presentation as well as for simulation purposes we prefer to write an single differential 
equation for the / r>n and another one for the / ri „, whose individual terms are switched on or off depending on the 
values of r, n. To this end, we have defined the following switch functions A, B,C, D, E and F: 



(96) 



dfr.n 

dt 



A(r, 


n) = 


1 


if 


r = 2kn = 


B{r, 


n) = 


1 


if 


r > 3&n > 1 


C(r, 


n) = 


1 


if 


(r = 2kn = 0)| 










(r > 3 & n < T(r) - 


D{r, 


n) = 


1 


if 


(r > 3&n < r(r- 1)) 


E(r, 


n) = 


1 


if 


(r = 2&n = 0)|| 










(3 < r < i?- 1) 


F(r, 


n) = 


1 


if 


r = 2&n = 



(97) 



(98) 

The above set of kinetic equations satisfy the constant density constraint which we monitor and verify in all the 
numerical simulations: 

L x + D 1 + Li + D 1 + 4P 

R r(r) 

+ ai + 2a 2 + ^ + 71 + 1 + fr,n) = const. 

r=2 n=0 

(99) 
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V. NUMERICAL RESULTS 



We arc interested in testing out the effectiveness of the proposed orientation and amplification cycle (Figure [l| as 
this is actually represented in the model. At the same time, we can also control the amounts of each type of additive 
and assess the relative importance of the two orienting effects, namely hydrophobic and kinetic. The results are 
quantified and presented in terms of a number of experimentally relevant chiral measures. The percent enantiomeric 
excess of the hydrophobic additives in solution is 

= {zjTlirT X 10 °- (1 ° 0) 
When hydrophilic amino acids are added, a corresponding percent enantiomeric excess can also be defined: 

ee « = mrm x m (101) 

We consider the percent excess of the occluded amino acids via 

occluded ^occluded ftr\n\ 

ee occ i ude (%) = — — x 100. (102) 

^occluded > ^occluded 

In keeping with the experimental reports, this counts the amino acids occluded in the host crystal enantiotopic faces 
exposed to the solution. An important aspect of the experiment is the degree of crystal orientation at the interface. 
The following parameter measures the orientation degree of the crystals at the interface: 

od{%) = t-—L x 100. (103) 

and gives us a direct measure of the orientation. Here / and / are the total amount of crystal of each type (i.e., 
pyramids and plates @) at the surface (omitting the floating inhibited seeds): 

R r(r) R F(r) 

/ = ^ fr,n f — ^ / r,n 

r— 3 n—0 r— 3 n—0 



The differential rate equations Eqs. (|87M95[) were numerically integrated using the version 7 Mathematica program 
package. The results were monitored to verify that the total system mass Eq. (I99|) remained constant in time. 
We organize and present our main results in terms of the types of additives employed, that is, hydrophobic and/or 
hydrophilic amino acids. 



A. Hydrophobic additives 

To simulate the following cases, only hydrophobic additives have been added to the saturated glycine solution. So 
we maintain zero concentration of the hydrophilic additives [ZJo = [Al]o = throughout the following sequence of 
simulations. All the numerical results have been obtained for the values R = 50, 7 = 0.1, a = 10~ 3 , 5 = 10 -6 , 
£ = 1, fio = 10~ 6 and f c — 0.01M. We provide some brief rationale for these selected values. The experiment reports 
that only a small fraction of the a- amino acids (0.02 — 0.2%) is subsequently occluded into the bulk of the growing 
glycine crystals. We therefore try to maximize the host crystal size R so as to be able to achieve small fractions of 
occluded guests. Computational limitations (memory and time) forced us to compromise and we thus choose R=50 and 
7 = 0.1 (for the desired values 7 = 0.0002 — 0.002, the corresponding R is so great as to exceed available computational 
resources). Next, we aimed to reproduce as closely as possible the reported details concerning the growth of crystals 
at the air/solution interface. The first step was to find reactions rates satisfying these conditions, as is shown in our 
Table 1. Since we do not have spatial information (and thus can not distinguish between pyramids and plates), we 
tried to find the reaction rates (at least the relation between them) reproducing this behaviour (supposing that pure 
platforms were the pyramids and the product of the inhibition yields the plates). This strategy led us to employ the 
values of a, S, £ and /i given above. 

We have used the same glycine concentrations as reported in the experiments*^. Thus, in a typical experiment 
10 g of glycine (m g i y = lOg) were dissolved in 30 ml (Vol = 0.03Z) of double distilled water, where its molar mass 
Mrrigiy = 75. 07g ■ mol -1 . Following this, the initial glycine concentration should therefore be [ajjo = m giy/(M' m giy x 
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Vol) = AAAmol ■ l^ 1 , and this is the value of initial glycine monomer concentration used in all the simulations 
reported here. For the glycine dimers we initially set [a 2 ]o = 0, and zero initial concentrations of the oriented glycine 
crystals at the interface: (f r . n )o = and {f r ,n)o = 0, respectively. Finally, from the initial enantiomeric excess 
eeo and the relative weight of additive to glycine, (w/w g i y ), we can deduce the initial concentrations [L]o and [Z?]o 
straightforwardly from the pair of Eqs. (|104[) . 

\L] 



(W/Wgly)m g ly(l + CCq) 

2Mm g i y x Vol 



rr)1 ( w /w g iy)m g iy(l - eep) ( 
[D] ° 2Mm glv x Vol ■ (104) 

We ignore the molar mass differences Mm between glycine and the ce-amino acid molar masses(Mm(Ala) = 
mmgmol- 1 , Mm(Leu) = Ul.Ugmol- 1 , Mm(Val) = 117.15gmol-\ Mm{Val) = 105.09.gmor 1 ). 



1. Racemic mixtures of hydrophobic additives 



. In this case, the composition of the complete solution is obtained by adding a racemic hydrophobic additive to 
the supersaturated solution of glycine. This means we start with zero initial excess of chiral hydrophobic monomers: 
ee = 0. We want to observe the effect of varying the relative mass of the chiral monomers with respect to the amount 
of glycine (w/w g i y ). The numerical results obtained in these simulations are shown in Table [IJ For each relative initial 
concentration of additives (w/w g i y ) we evaluate the total concentration / + / of the oriented mixed host plus guest 
crystals at the interface summed over both orientations, and P, the total concentration of nuclei or seed crystals killed 
by the kinetic inhibition effect and their associated chiral polarizations or excesses. Recall that both the hydrophobic 
and kinetic effects are operative for hydrophobic amino acid additives. 

As we see in TabkUl for all cases the interface is covered with crystals of both orientations. This is not surprising, 
since we start from a racemic composition, but no initial "Adam" crystal, this must lead to an equal ratio / : / = 1 : 1 
of the enantiotopic faces exposed to the solution, thus od{%) = 0. By the same token, the ee of the solution and 
the net ee oc of occluded monomers must also be zero. There is no symmetry breaking in this idealized situation. 
In the actual experiments, floating glycine crystals can exhibit two distinct morphologies, either pyramidal or plates 

||. The former are the result of the hydrophobic orienting effect, while the latter morphology results from the 
kinetic inhibition effect. Because of this, we can infer the crystal morphologies implicit in our simulations: indeed, 
the concentrations / and / correspond to the pyramidal form (resulting from the hydrophobic effect) whereas the 
concentration P corresponds to the plate-like form (resulting from the kinetic effect). 

Thus for values of w/w g i y as low as 10~ 4 , we found enantiomorphous pyramids showing thus no macroscopic 
evidence for the kinetic retardation of growth by for example, leucine. By increasing the additive concentration, 
a morphological change of the floating glycine crystals is numerically observed, from pyramids / + / to plates P 
(inhibited seeds). The concentration of the pyramids decreases while that of the plates increases as we increase the 
concentration of the racemic additive. This simulation result is in qualitative accord with the experiment where an 
evident morphological change from pyramids to plates was observed upon increasing the leucine concentration (i.e., 
the hydrophobic additive employed there). @. The initially racemic amino acid solution remains racemic at all later 
times, no matter how large the initial amount of additives is. 



TABLE I: Orientation of floating glycine crystals in the presence of mixtures of (L,D) hydrophobic additives 



Hydrophobic 


/ + / 






P 






od 


ee 




W/lVgly 


mol ■ 


r 1 




mol ■ 


r 1 




(%) 


(%) 


(%) 


1CT 4 


8.15 ■ 


10" 


-5 


4.38 


10" 


-6 











icn 3 


6.33 • 


10" 


-4 


3.54 


10" 


-4 











10 -2 


2.38 • 


10" 


-3 


1.51 


10" 


-2 











lO" 1 


3.48 • 


10" 


-3 


2.12 


10" 


-1 












2. Racemic mixture hydrophobic additives + an initial oriented glycine crystal 
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The reaction scheme in Figure [JJ starts from a first random crystallization, so we include this feature in the initial 
conditions. We consider the previous situation (ee = 0) and suppose that there is one single crystal of glycine exposing 
say the (010) face toward the solution at an initial time; this means starting with a small initial concentration of 
(/2,o)o- Since only one initial crystal can not be taken, it can however be "modeled" taking by a fraction of the 
critical interface concentration f c of glycine crystals: (/2,o)o = x f c , where x « 1 is given in the first column of 
Table [Til There we display the numerical results for this particular situation, starting from an initial concentration 
of (LD)-hydrophobic amino acids w/w g i y = 0.01. The simulation results show a preferential (010) orientation of the 
glycine crystals at the surface induced by the initial (010) crystal. Even if this preferential orientation is not exclusive 
od < 100%, an optically active solution is still generated, ee = 100%. 

The single crystal occludes the D enantiomers enantioselectively, thus enriching the aqueous solution with the L 
monomers. If this small excess can induce preferential orientation of the further growing crystals of glycine again 
with the (010) face pointing toward the solution, replication will ensue by cascade mechanism finally leading to a 
separation of enantiomeric territories [6]. From Tabic UJ] we see that increasing the initial concentrations of the first 
oriented Adam crystals leads to increased degree of orientation od, and to increased enantiomeric excesses of both the 
amino acids in solution ee and of those that have been occluded by the pyramidal crystals, ee oc - Dynamical aspects 
of the symmetry breaking can be appreciated from the time evolution of the concentrations and chiral excesses as 
displayed in Figure [3] 

TABLE II: Orientation of floating glycine crystals in the presence of mixtures of (L,D) hydrophobic additives and "one" initial 
oriented crystal 



(/a,o)o//c 


/ + / 


P 




od 


ee 


GCoc 


X 


mol ■ i" 1 


mol ■ 


r 1 


(%) 


(%) 


(%) 


0.01 


3.02 • 10" 3 


1.51 


10" 2 


30.04 


100 


7.10 


0.05 


4.59 • 10~ 3 


1.47 


10 -2 


62.46 


100 


21.55 


0.1 


6.11 • 10" 3 


1.42 


10" 2 


76.32 


100 


33.60 



100 



^ 80 




time 

FIG. 3: Racemic mixture hydrophobic additives + an initial oriented glycine crystal. Initial conditions: eeo = 0, w/w g i y = 0.01, 
[/2,o]o = 10 -3 f c , ee h — 0, and Wh/w g i v = 0. Upper figure: dotted line represents the enantiomeric excess ee of solution and the 
solid line the orientation degree od. Lower figure: the dotted lines represent the concentrations of the hydrophobic enantiomers 
(where [L] is the upper dotted curve), the solid lines represent the total concentration of oriented crystals (where [f] is the 
upper solid curve). 
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3. Chiral hydrophobic additives 



Instead of racemic hydrophobic additives we next consider adding a chiral hydrophobic compound (i.e., L-a-amino 
acids) to the supersaturated solution of glycine. This means we now start with an initial excess of the chiral hy- 
drophobic monomers in solution so that eeo = 1. Table Hill shows the results for different concentrations in terms of 
the quoted w/w g i y values. As expected, in the presence of a resolved hydrophobic a-amino acid, the system achieves 
exclusive orientation at the interface: od — 100% and the excess ee of the remaining amino acids in solution is 100%. 
In this situation, there can be no amino acids occluded by the oriented crystals at the interface since the system 
lacks the Z?-monomer, hence ee oc is not defined. Moreover without D monomers, the kinetic effect is inoperative and 
so P = 0: no plates can be formed, only pyramidal crystals. When L-a-amino acids are used, these pyramids are 
exclusively (010) oriented (the face exposed towards the solution). By symmetry, D amino acids would induce the 
enantiomorphous (OlO) oriented pyramids^- 



TABLE III: Orientation of floating glycine crystals induced by resolved L-hydrophobic additives (See Table 111)0]) 



(L) 

Hydrophobic 

W/Wgly 


f + f 
mol ■ I -1 


P 

mol ■ I' 1 


od 

(%) 


£6 &&oc 

(%) (%) 


0.01 


10" 2 





100 


100 


0.05 


10~ 2 





100 


100 


0.1 


1.01 • 10" 2 





100 


100 


0.5 


1.02 • 10" 2 





100 


100 



B. Hydrophilic additives 



We now treat the general situation where both hydrophobic and hydrophilic additives are added to the super- 
saturated glycine solution. The following numerical results have been obtained for the values R = 50, 7 = 0.1, 
a = 10~ 3 , 6 = 10~ 6 , /3 = £ = 1, fxo = 10 -6 and f c = 0.01M. As before, the initial glycine monomer concentration is 
[ai]o = 4.44 moZ I and for the initial glycine dimer concentration we set [aa]o = and zero initial concentrations of 
the oriented glycine crystals at the interface: ( f r , n )o = 0, (/ r ,n)o = 0. Controlling the amount of hydrophilic additives 
gives us an independent control of the kinetic inhibition effect. 



1. Racemic mixtures of hydrophilic additives 



The complete amino acid solution is now composed of a mixture of both racemic hydrophobic a-amino acids 
and by racemic hydrophilic a-amino acids. Here, the role of the hydrophobic additive is only to induce oriented 
crystallization at the air/water interface. Since its composition is racemic, this leads to an equal ratio / : / = 1 : 1 
of enantiotopic faces exposed to the solution, and for this we use a concentration corresponding to w/w g i y = 0.01. 
Thus the corresponding initial excesses of chiral hydrophobic and hydrophilic monomers are eeo = and ee^ = 0, 
respectively. There is no initial Adam crystal. We then observe what effect, if any, varying the amount of the 
hydrophilic additives has on this initial orientation. 

Numerical results showing the effect of varying the relative mass of hydrophilic chiral monomers (w/w g i y ) are 
shown in Table IIVI No matter how large the amount of hydrophilic monomers, in the presence of racemic additives 
of both types, the surface is still covered with crystals of both orientations which immediately implies that the 
net excess of occluded monomers is zero. Note the solution achieves an extremely feeble optical activity. In fact, 
these excesses are only slightly greater or at the same level as the enantiomeric excess expected from pure statistical 
fluctuations 13[ • Just as reported in the experiment, we observe how the amount of plates is much greater than the 
amount of enantiomorphous pyramids: P > (f + /)[6|. 
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TABLE IV: Orientation of floating glycine crystals in the presence of racemic mixtures of hydrophilic additives and 1% 
hydrophobic (L,D) amino acids 



Hydrophilic 
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2. Chiral hydrophilic additive 



Just as in the previous simulation, the amino acid solution is composed of a racemic mixture of hydrophobic a-amino 
acids merely to induce the X : Y = 1 : 1 crystallization at the air/water interface and again we use w/w g i y = 0.01. 
But this time, we add chiral hydrophilic additives, and without loss of generality we take the L enantiomer. This 
is then a situation described by the initial excesses of hydrophobic and hydrophilic monomers given by eeo = and 
eeho = 1, respectively. Here, we are interested in the effect of varying the relative mass of hydrophilic chiral monomers 
(w/uigiy). Table IV1 presents the numerical results for this situation. As we can see there, the presence of a chiral 
hydrophilic additive induces a preferential orientation, it is an indirect mechanism, since the kinetic effect inhibits 
the growth of the crystal nuclei that would otherwise expose the V-face towards solution. Even if the orientation is 
not exclusive but only preferential, the solution does become optically active for low concentrations of the hydrophilic 
additive. And from the Table we see that territorial separation of the enantiomers is achieved. 

TABLE V: Orientation by kinetic effect of floating glycine crystals in the presence of resolved hydrophilic additives and 1% 
hydrophobic (L,D) 
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Hydrophilic 
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P 






od 


ee 




W/Wgly 


mol ■ I -1 




mol 


r 1 




(%) 


(%) 


(%) 





2.38 • 10" 


:i 


1.51 


10" 


-2 











10" 6 


2.38 • 10" 


3 


1.51 


10" 


-2 


0.51 


100 


0.08 


10 -5 


2.51 • 10" 


3 


1.51 


10" 


-2 


8.43 


100 


1.56 


10" 4 


3.28 • 10" 


-3 


1.52 


10" 


-2 


38.88 


100 


9.21 


10" 3 


6.10 ■ 10" 


-3 


1.52 


10" 


-2 


79.03 


100 


33.36 


10~ 2 


1.00 ■ 10" 


-2 


1.51 


10" 


-2 


96.26 


100 


76.09 


lO" 1 


1.00 • 10" 


-2 


1.51 


10" 


-2 


99.51 


100 


96.60 



The hydrophilic additives have an indirect orienting effect upon the floating glycine crystals in that they inhibit 
or kill the glycine nuclei so that these are unable to occlude. Only the hydrophobic additives can directly induce 
orientation. A typical time evolution of the symmetry breaking in this situation is depicted in Figure 2] 

3. Hydrophobic vs. Hydrophilic additives 



Next we study how hydrophobic and hydrophilic additives of opposite chiralities compete to induce a preferential 
orientation of the glycine crystals at the interface. To emphasize this competition, the amino acid solution will be 
composed of strictly (L)-hydrophobic and (D)-hydrophilic a-amino acids. So, initially we have eeo = 1> eeho = — 1- 
Take the (L) hydrophobic additive concentration to be w/w g i v = 0.01, and we vary the (D) hydrophilic additive 
concentration. Table IVII shows there is a clear ability of the resolved hydrophobic a-amino acids to induce a specific 
orientation of the floating glycine crystals, even in the presence of large excesses of the hydrophilic a-amino acids of 
the opposite absolute configuration p. This can be compared with Table IV 6], which established experimentally the 
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FIG. 4: Chiral hydrophilic additive. Initial conditions: eeo = 0, w/w g i y = 0.01, eenO = 1, Wh/wgiy — 0.01. Upper figure: the 
dotted line represents the enantiomeric excess ee of the hydrophobic additives in solution, the solid line the orientation degree 
od, and the dashed line is the eeu of the hydrophilic additive. Lower figure: dotted lines represent the concentrations of the 
hydrophobic additives (where [L] is the upper dotted curve), the solid lines represent the total amount of oriented crystals of 
each type ([f] is the upper solid curve), and dashed line represents the hydrophilic monomer concentration: [L]. 



dominance of the hydrophobic effect over the kinetic effect. Note: in our model, ee oc is zero in this situation since 
there are no D hydrophobic monomers that would otherwise be occluded. 



TABLE VI: Orientation of floating glycine crystals by hydrophobic effect in the presence of resolved hydrophilic additives and 
1% (w/wgiy) hydrophobic amino acids of the opposite configuration 
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4- The amplification step 



The hydrophobic resolved a-amino acids are advantageous for such a study with respect to the hydrophilic ones 
because both the kinetic inhibitory effect and the stabilization of the nuclei by hydrophobic effect act in the same 
direction for orientation of the growing floating crystals of glycine. For this reason, we next study the crystallization 
of glycine in the presence of partially enriched mixtures of (L,D)-hydrophobic amino acids at various concentrations. 

The presence of an excess of hydrophobic (D)-amino acids will favor the (010) oriented nucleation at the interface 
while at the same time preventing the growth of (OlO) oriented crystals from the solution. The experimental [f| 
correlation between the initial enantiomeric excess of the solution with the total concentration needed to obtain 
complete orientation of the floating glycine crystals is reproduced here in Fig|5]for reference. 
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The correlation between the initial enantiomeric excess of the solution with the total concentration needed to obtain 
maximum orientation of the floating glycine crystals is shown in Fig. [SJ In the Table I VIII we list all the od values 
obtained for the different initial enantiomeric excesses and concentrations, the numbers in boldface correspond to the 
maximum orientations and these are plotted in Fig. [5] The point we wish to illustrate here is simply that our model 
succeeds in capturing the general trend observed in the experiment, namely that smaller initial enatiomeric excesses 
require greater initial hydrophobic amino acid concentrations in order to achieve a (maximal) crystal orientation. 
Figure [7] shows the temporal evolution of this situation. 

100 I — i — i — i — i — . — . — i — i — i — i — . — i — . — i — i — i — . — . — . — i — i — I 
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FIG. 5: Correlation between the initial leucine enantiomeric excess in solution and the total concentration needed for the 
complete (010) orientation of the floating glycine crystals. From reference [g]. 
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FIG. 6: Correlation between the initial enantiomeric excess of the hydrophobic additive in solution and the total concentration 
needed for the maximal (OlO) orientation of the floating glycine crystals. See Table I VIII and text for discussion. 

Experimental proof that hydrophobic effect plays a dominant role in the orientation of the glycine crystals can 
be deduced from Table IV of reference @. Nevertheless, both effects are always present for the hydrophobic amino 
acids: those in solution contribute to the inhibitory kinetic effect, this is depicted in Figure [5j This also explains why, 
at higher concentrations, the glycine plates observed in the experiment are so thin. At lower initial ee, the higher 
concentrations needed appear to forbid formation of floating plates altogether Q. 

VI. CONCLUSIONS AND DISCUSSION 

We have presented a mathematical model for achieving resolution of racemic solutions of alpha-amino acids and 
glycine into enantiomeric territories based on a chemical scheme proposed by the Rehovot group some years ago. Their 
crystallization experiments provide a simple model for the generation and amplification of optically active amino acids 
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TABLE VII: Orientation degree od of floating glycine crystals in the presence of enriched mixture of hydrophobic additives. 
For each initial enantiomeric excess, a maximum degree of orientation is reached, to identify it we first locate the pair of values 
(the range) containing this maximum (bold font) 
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TABLE VIII: Maximum orientation degree od of floating glycine crystals in the presence of enriched mixture of hydrophobic 
additives for each initial enantiomeric excess 
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in prebiotic conditions. The glycine and alpha amino acid system may be relevant to the origin of optical activity 
as it involves compounds that are among the simplest building blocks of life. In this vein it is interesting to point 
out that a recent prebiotic synthesis of protobiopolymcrs under alkaline ocean conditions has yielded a variety of 
amino acids with glycine being the most abundant [9]. In an astrophysical context, amino acids have been detected 
in room-temperature residues of UV-irradiated interstellar ice analogues, and again glycine was found to be the most 
abundant amino acid[l6j. 

Centrosymmetric glycine crystals were employed as substrates for the total separation of occluded alpha-amino acids 
into enantiomeric territories. These amino acid additives are occluded enantioselectively through the enantiotopic 
(010) and (010) faces of the glycine crystals. Such crystals when floating at the air/solution interface, and if properly 
oriented, may incorporate only one of the two enantiomer additives present in the solution. Complete (010) or (010) 
orientation is induced by both a kinetic and a hydrophobic effect. The former effect acts through an inhibition of the 
nucleation and growth of that enantiomorph which interacts with the resolved additive. The latter effect is due to 
the induction of a specific enantiotopic face orientation (e.g., (010) exposed toward the solution) by the hydrophobic 
amino acids at the interface. Combination of both effects acts in the same direction driving exclusive orientation of the 
glycine crystals and thus triggers an amplification starting from an initial random oriented crystal and a solution with 
low initial enantiomeric excess. This is represented schematically by Figure [TJ Our mathematical model results from 
translating this scheme into basic reaction steps and then into a corresponding system of differential rate equations 
which we then simulate to test out the hypothesized mechanisms and to underscore the salient features of the original 
experiment. 

Our immediate goal here is to capture the essential mechanisms responsible for the symmetry breaking in the sim- 
plest terms possible. This has led us to make certain simplifications, in keeping with this aim. A major simplification 
is to adopt reaction rates that are independent of the instantaneous size of the glycine crystal and of the number of 
occluded amino acid guest monomers. Another simplification is to model the processes without using explicit coordi- 
nate dependence, which otherwise would have necessitated the introduction of partial derivatives, diffusion constants, 
spatially dependent concentrations and a coordinate system that distinguishes the bulk three dimensional solution 
from the bounding two-dimensional air/water interface or layer. There is no doubt that these spatial aspects play 
a supporting role in the actual experiment, but they are not the primary cause of the mirror symmetry breaking 
observed there. The fact that the coordinate-free model leads to chiral symmetry breaking and to the territorial 
separation of the enantiomers is proof of this. We have provided an analytic linear stability analysis which gives inde- 
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FIG. 7: Partially enriched mixtures of (L,D)-hydrophobic amino acids. Initial conditions: eeo = 0.1, w/w g i y — 0.01, ee^O = 0, 
Wh/w g iy — 0. Upper figure: the orientation degree od (solid line) and the enantiomeric excess ee in solution (dotted line). Lower 
figure : concentration of the hydrophobic monomers (the upper dotted curve is [L]). Solid lines represent the concentration of 
oriented crystals ([f] is the upper curve). 



pendent confirmation of this. The mathematical model leads to separation of enantiomeric territories, the generation 
and amplification of optical activity by enantioselective occlusion of chiral additives through chiral surfaces of glycine 
crystals. 

We emphasized earlier the experimentally remarkable feature that the chemical system does away with the need for 
mechanical energy for achieving the amplification of the ee of the amino acids in solution. To this end we can identify 
an aspect of the experiment that serves as an effective " driving force" that maintains the system out of equilibrium, a 
seemingly necessary condition for the (permanent |24j |) breaking of mirror symmetry. Indeed, the experiment requires 
supersaturated solutions of glycine. To achieve supersaturation, the system must be (i) cooled down and/or (ii) the 
solvent must evaporate. This then adds the element of irreversibility. 

Recent works related to the original experimental model we consider here are summarized in a current review 
of the role of crystalline architectures as templates relevant for the origins of homochirality [12]. In brief, similar 
stochastic separation of a-amino acids has also been achieved when racemic alpha-amino acids are occluded within 
enantiomorphous (3- glycine crystals [l7j . When this form is grown in either porous materials [l8l [l9j or small solution 
volumes[2(|, the beta polymorph glycine crystallizes as long needles. When the /3-form is grown in the presence of 
DL-amino acids, the L-molecules are occluded only in one of the /3-glycine enantiomorphs, while the D's are occluded 
only in the other. Thus, if a small number of one the j3 enantiomorphs crystallizes first, they will occlude only one of 
the alpha amino acid enantiomers, thereby enriching the solution in the other amino acid. This excess should prevent 
the nucleation of fresh /3-Gly crystals of the opposite handedness and thus lead to an enrichment of the ee of the 
amino acids in solution. 

The present model is much simplified with respect to the actual experimental situation in solution. It is therefore 
worthwhile to identify potential future efforts aimed at "dressing" this minimal model so as to increase its relative 
"complexity" , only of course when this is warranted or justified experimentally. One such effort we are currently 
investigating is a more realistic modeling of the self-assembly of the hydrophobic amino acid monolayer at the 
air/water interface and its role in precipitating the oriented glycine crystals. Further theoretical work aimed at 
incorporating other physical and chemical aspects to the underlying scheme could also include, for example, mass 
transport and diffusion, and a detailed modeling of the the glycine crust formed at the interface. We hope to report 
results along these lines elsewhere. 
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